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Abstract 

Quantum theory for measurements of energy is introduced and its conse- 
quences for the average position of monitored dynamical systems are ana- 
lyzed. It turns out that energy measurements lead to a localization of the 
expectation values of other observables. This is manifested, in the case of 
position, as a damping of the motion without classical analogue. Quantum 
damping of position for an atom bouncing on a reflecting surface in pres- 
ence of a homogeneous gravitational field is dealt in detail and the connection 
with an experiment already performed in the classical regime is studied. We 
show that quantum damping is testable provided that the same measurement 
strength obtained in the experimental verification of the quantum Zeno effect 
in atomic spectroscopy [W. M. Itano et al., Phys. Rev. A 41, 2295 (1990)] is 

made available. 
03.65.Bz, 42.50.Wm 
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I. INTRODUCTION 



A number of new experimental techniques involving devices with noise figures close to the 
minimum dictated by the uncertainty principle has originated the demand for quantitative 
predictions of quantum measurement theory fT]]^]. One of the most important steps along 
this path was made by Itano et al. ||: they experimentally showed that the effect of the 
measurement on the observed system leads to a freezing of its free dynamical evolution, 
the so-called quantum Zeno effect In such a case one is looking at the occupancy 

probability of one level, that happens to depend upon the process of measurement itself. 
More general quantum measurement effects can be obtained by measuring one observable and 
by looking at the dynamical evolution of another observable, as it has been already discussed 
for a mesoscopic structure in |J. This agrees with the original spirit with which Sudarshan 
opened the debate on the quantum Zeno effect and according to which measurement of the 
position of a decaying particle should influence its lifetime, i.e. another observable quantity. 

In this paper we discuss a model for a system subjected to quantum measurements of 
energy including their influence on other observables, specifically position measurements. 
The act of measurement affects the average position of the system giving rise to an average 
localization of the motion, a quantum damping without classical counterpart. The paper is 
organized as follows. In Section II we introduce the general formalism for measurements of 
energy distinguishing between nonselective and selective cases and we evaluate the average 
position of the measured system in both schemes. The example of the harmonic oscilla- 
tor, for which analytical evaluations are possible, is also dealt. In Section III a particle 
in a homogeneous gravitational field bouncing on a perfectly reflecting surface is analyzed 
in detail. In particular, the dependence of the quantum localization effect upon the rele- 
vant parameters is studied both by means of numerical tools as well as in the semiclassical 
limit, by exploiting the WKB approximation. In Section IV possible scenarios to look for 
the predicted quantum damping are discussed with particular emphasis on the previously 
considered particle bouncing on a reflecting surface. We focus on a case very close to the 
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experimental situation already reported in [[/J, giving a discussion of the parameters that 
have to be reached to look for the predicted effect. 



II. QUANTUM MEASUREMENTS OF ENERGY: GENERAL FORMALISM 

Since the original proposal of Von Neumann || , important progresses have been made 
in understanding quantum measurement theory. In particular, it has been recognized that 
the dynamics of the system in presence of the measurement undergoes modifications, with 
respect to the closed system dynamics, that can be taken into account by means of an ef- 
fective master equation or a Schrodinger equation for mixed and pure states respectively. 
Thus the original doubling of dynamics, a free Schrodinger evolution during the non mea- 
surement periods and an abrupt state collapse during an ideal, instantaneous measurement, 
has been replaced by a unique dynamical approach. In this context, some models represent 
the meter as a particular environment that interacts with the observed system and extracts 
information from it. The dynamics of a system interacting with an environment is most con- 
veniently described in terms of a reduced density matrix operator pit) , obtained by tracing 
out the degrees of freedom of the environment from the density matrix operator of the entire 
system (marginalization procedure). The unitary evolution of p(t) for the isolated system is 
modified to an irreversible one due to the interaction with the environment. In the limit of a 
Markovian environment, a dynamical law described by a completely positive semigroup can 
be postulated for the evolution of the open system PJlQfl, resulting in the following master 

equation for the reduced density matrix operator: 
d i 1 n 

Jt p(t) = -j- [H(t),p(t)} +-E ([Lv(t)p(t),Lv(t) 1 } + [Ut),p(t)L u (ty]) (1) 

where H (t) = H(p, q, t) is the Hamilton operator for a general nonautonomous system and 
L u (t), v = 1, ... ,n, are the so-called Lindblad operators, that are supposed to model the 
effects of the environment on the system. The above equation describes the general case of a 
quantum open system. The evolution of a quantum system subjected to a measurement pro- 
cess corresponds to the particular case where the environment is the measurement apparatus 



and the Lindblad operators are proportional to the measured quantities. If the measure- 
ment of a single observable represented by the operator A(t) = A(p, q, t) is considered, the 



corresponding Hermitian Lindblad operator can be chosen to be L(t) = L(t) ] = ^K(t)A(t). 
The function n(t) has dimensions [k] = [t _1 A^ 2 ] and represents the coupling, in general 
time-dependent, of the monitored system to the measurement apparatus. The result of the 
measurement is 

<t) = Tr (A(t)p(t)) > (2) 

where overlining denotes a statistical average over individual results associated to pure states 
that form the incoherent mixture described by the reduced density matrix operator p(t). 
By moving to the coordinate representation of the reduced density matrix operator, 

p(quQ2,t) = p(t) \q 2 ) , (3) 

the evolution equation obtained from ([I]) is 



A ( ~^^~~> 9i> t) - A q 2 , t 



\p{qx,q2,t). (4) 



Equation @ gives the general description of a system in which the observable A is 
continuously monitored and the result of the measurement is not known in advance, the 
so-called nonselective measurement process. 

If continuous measurements of energy characterized by a constant coupling are con- 
sidered (A = H), Eq. (4) for the evolution of the density matrix specializes to 

Pnm ~j~i.E n -E771) 'Pnm ^ ("^1 ^m) Pnm 1 (5) 

where 

Pnm{t)= / dq dq' p(q,q',t)(j)* n (q)(j) m (q') (6) 



and {4> n } is the complete set of energy eigenstates, corresponding to the energy eigenvalues 
E n . It is easy to show that the solution of the equation (H) can be written as 



Pnmit) = exp ^~{E n - E m )t -^-(E n - E m ft^ p nm (0) , (7) 
and consequently 

p(q,q',t)=J2pnm(t)<i> m (<l')M<l)- (8) 

nm 

Notice that the effect of the measurement in ([?]) is to diagonalize the density matrix after 



enough time, the so-called decoherence induced by the measurement ||11|| . 

Without loss of generality, we can restrict our analysis to the position operator Q. The 
average position at time t is then evaluated as 

(Q(t)> = Jdq dq' p(q,q',t) (q'\Q\q) = J dq p(q,q,t) q (9) 

and therefore from (|8]) 

(Q(t)) = T,Pnm{t)(n\Q\m) . (10) 

nm 

By inserting (|7|), the average position in the nonselective case we have just considered is 
finally written as 

(Q(t)} = E ex P {~U E n ~ E m)t -^(E n - E m ) 2 t\ P nm(0)(n\Q\m) , (11) 

nm ill ^ ) 

where the effect of the measurement leads to exponential decaying behaviour with state- 
dependent time constants 
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It is worth to observe that the off-diagonal terms of the sum contributing to the average 
position vanish in the asymptotic limit t — > oo when > 0. As a consequence, the 
average position tends to be localized around the stationary value due to the diagonal 
contributions. From the conceptual viewpoint, this effect is nothing but the manifestation, 
in the configuration space, of the previously mentioned decoherentization process. 

Another interesting situation arises when the energy measurement is already performed 
and the result is known. The theory should be able, in this case, to complete the knowledge of 



the system by evaluating the corresponding wavefunction conditioned to the known result of 
the energy measurement, the so-called a posteriori selective measurement. In 0] it is shown 
that, starting from (|j), an effective Schrodinger equation for an a posteriori measurement is 
obtained for the restricted wavefunction ip[ a ] (q, t) as: 



ih 



dlp[a] (g, t) 

dt 



d 

H(-ih—, q, t) - 1% K[ a ] (t) 



d 

A(-ih—,q,t) - a(t) 



1 2 



ip [a] (q,t) , (13) 



where we have denoted by a the result of the measurement and by [a] the associated func- 
tional dependence. Equation (|13|) formally holds also when the result a(t) is unknown and 
has to be predicted with probability distribution ||^[ a ](t)|| 2 (a priori selective measurement); 
a different interpretation of (|T3|) is however required, for which we refer to fl2fl . Coming 



back to continuous measurements of energy with a constant result E(t) = E, Eq. ([13]) has 
the solution H|l3"| 



4>[E\(q,t) = E c «(°) ex P {-iE n t/h - n E (E n - Eft) <j> n (q) . (14) 

n 

Here, as before, (f> n (q) are energy eigenfunctions, c n (0) the corresponding initial projection 
coefficients, and normalization to unity is lost due to the non Hermitian character of Eq. 
([1"3|), that expresses in more physical terms the branching of the wavefunction among all the 
possible alternatives. By exploiting (fH|), the average position in this selective case, hereafter 
denoted by (Q{t))[E\, is therefore written as 

( Q(t)) [E] = {M\Q\M)I(M\M) = [E |c fc (o)| 2 ex P {-2^( J E; fc - Eft}] 1 ■ 

L k J 

• J2c n (0)cl(0)exp{-i(E n -E m )t/h-K E [(E n -E) 2 + (E m -E) 2 ]t} (m\Q\n) . (15) 

nm 

The link between the average position in the nonselective (0) and selective (|i~5|) cases is 
established, in the usual way jl2| , by summing over all the possible selective measurement 
processes, i.e. 

(16) 



P(q, q', t) = J d[E] ij)* [E] (q, t) ip [E] (q 1 , t) , 

that also shows how the nonselective measurement, described in terms of the density ma- 
trix, can be seen as a functional integration over all the possible selective measurements, 
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represented by the restricted wavefunction. Note also that the decay time constants in ( |i"5|) 
explicitly depend upon the value E registered during the measurement. 

The theoretical scheme we have presented can be easily implemented for a simple system 
as the harmonic oscillator. Let us denote by m and u the mass and the angular frequency 
respectively. The effect of the measurement on the average position in the nonselective case 
can be guessed by simple inspection of Eq. (|TT|). Indeed, due to the equal spacing between 
adiacent levels, one immediately realizes that the measurement damping effect factors out 
and acts on the remaining unmeasured evolution as a purely exponential damping with an 
unique time constant r = 2/KEh 2 LU 2 , that is: 



(Q(t)) Ki 



exp 



(Q(t)) 



K E =Q 



(17) 



In the selective case instead the corresponding expression fll4|) can be explicitly computed 

if, as usual, creation and annihilation operators a, are introduced so that 

/ h 

(18) 



Q 



2muj 



a + a 



Then, owing to the well known selection rule for the matrix elements of the position operator 
between energy eigenstates, namely 

/ h 

(19) 



(n\Q\m) 



h 



(n\(a + a})\m) 



j 2mio 1 V 2mw 

the following expression for the average position is found: 

' h \ 1/2 



Vm8 n , m -i + Vm + 15 n , m+ i) , 



«?(*)> 



2muj 



s -1 oo 

exp[— 2K E {E k — E) 2 t] \ ■ ^ Vn + 1\ (20) 

t > cos(ut + 8) , 



\Ck 



■ exp < — 2ke 



1 



n 



i 2 



hu-E 



+ hu 



n H — I Yiuj — E 
2 



+ 



n=0 



where 9 is the relative phase between c n and c n+ \. Unlike the nonselective case, no factor- 
ization of the damping factor is allowed here. By referring to the decay constant for the 
nonselective measurement r, we get 



T\E](n 



T 

4 



1 E ) I E 

2 Tiuj J hto 



n -1 



(21) 



The time constants will now depend upon the registered energy E, decreasing with the 
difference between the ratio E/hu and the average number of quanta in the state. 
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III. FREE FALL OF A PARTICLE AND BOUNCING ON A REFLECTING 



SURFACE 

Another stimulating example is the monitoring of the energy of a particle falling in an 
homogeneous gravitational field and bouncing on an elastically reflecting surface. Experi- 
ments already performed have shown that multiple bouncing of atoms on surfaces is possible 
J7[ and the classical sources of damping have been understood ffM . The experiment has been 



performed by measuring the average of an ensemble of independent atoms and therefore the 
nonselective approach described in Section II is more adequate. 

Let us consider the height z over the surface and the potential energy 

{mq z z > , 
(22) 
+00 z < . 

The Schrodinger equation for positive z and energy eigenvalue E is 

-£S + <^-^ = < 23 > 

and we assume, consistently with the boundary conditions of the problem, ip(z) =0 for 
z < 0. Equation (^) is exactly solvable and the solutions are expressed in terms of the Airy 
functions [|T^ 

<f> n (z) = C n Ai {z/zq - A n ) , (24) 
where we have introduced the characteristic length of the system 



/ tr \ m 



The CnS are normalization constants and — A n is the n-th zero of the Airy function, related 
to the corresponding energy eigenvalue through the following relation: 

According to these conventions, z X n represents the n-th classical turning point, also denoted 
by z n . 



The distance between consecutive zeroes of the Airy functions |A n — A n+ i| — > forn — > oo. 
It is therefore recognizable that, for this system, unlike the harmonic oscillator, the quan- 
tum damping is large for states formed by low-energy pairs of consecutive eigenstates and 
vanishes in the (classical) limit of high energy. To understand the mechanism of damping we 
have numerically studied various cases corresponding to different initial preparations of the 
system. The numerical accuracy of the program, tested with the already analytically solved 
harmonic oscillator, is of the order of 0.1% for reasonable values of the space-time lattice. 
If energy eigenstates are considered, the average position is constant in time but, unlike the 
case of the harmonic oscillator, it is different from zero due to the nonvanishing diagonal 
matrix elements of the position operator. However, these stationary states are not affected 
by the presence of a measurement coupling. On the other hand, when a superposition of 
two energy eigenstates is assigned at zero time, the average position will harmonically oscil- 
lates between a minimum and a maximum value with angular frequency uo = (E n — E m )/h 
(n = 2, m = 1 in the example of Fig. la). By comparison, the constant values of the average 
position in the two energy eigenstates are also shown, the quantities 

Qn = J dzMz)* z Mz) ! Q22 = J dzfa(z)*z<h(z) , (27) 

which allows one to write the average position of their superposition ip = c\4>i + c 2 4>2, as 

(Q(t)) = |ci| 2 Q n + |c 2 | 2 Q22 + 2|ci||c 2 |<3i2 cos 

with 

Q 12 = Jdz<p 2 (z)* z faiz) = Q* 21 . (29) 

In presence of a measurement a pure exponential damping arises that still reminds the 
oscillatory behaviour (as shown in Fig. lb) or that is overdamped (as in Fig. lc), depending 
upon the measurement coupling through the time constant T\i = 2 [ke{E^ — -Ei) 2 ]" 1 . The 
time development is centered around the mean value of the unmeasured evolution ( f28|) . The 
effect of the energy measurement vanishes in the limit of large quantum numbers, namely 
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(Ei - E 2 )t 

h 



+ 9 



(28) 



the classical limit. This can be easily shown by exploiting the WKB approximation, that 
gives the spacing between two consecutive levels as reported in fHf 

E n = ^(27) V2 • (30) 

Since it is difficult to obtain superpositions of only two energy eigenstates, the attention 
must be focused on states that are more likely to be produced when an atomic cloud is 
prepared. This is the case of the Gaussian states, already extensively studied for being the 
quantum states closest to classicality, although still far from representing an atomic cloud. 
As discussed in [O] , Cesium atoms in magnetooptical traps lead at best to an initial radius 
of the cloud equal to z,- t ~ 50 fim with an rms momentum corresponding to a velocity spread 
of ~ 2cm/s, leading to an uncertainty product Zip Zi ~ 2000 h. Nevertheless, Gaussian 
states constitute entities simple enough to be analyzed, allowing a physical understanding 
of their behaviour, and at the same time complex enough to maintain the typical features 
of the more realistic situations. Indeed, the striking difference from the already analyzed 
case is that now the energy eigenstates expansion cointains different eigenstates and a more 
complex dynamics is obtained also in the unmeasured case (Fig. 2a). When the effect 
of the measurement is taken into account, a damping of the motion is obtained in all the 
explored cases, although it is not a pure exponential damping, as shown from the weak 
reformation of the damped oscillations in Fig. 2b. By varying the height of the atomic 
center of mass for a Gaussian state the relevant contributing eigenstates also change. In 
general, as also intuitively understandable, for a height z the greatest contributions will 
stem from the eigenstates having the classical turning point zo\ n closest to z. Moreover, the 
number of contributing eigenstates will depend upon the width of the initial Gaussian state. 
In Fig. 3 a configuration similar to Fig. 2 but with a larger width of the Gaussian state 
is depicted. In this case, more energy eigenstates significantly contribute to the expansion 
and in particular the lower ones determine, due to their larger energy difference, a faster 
dynamics for the average position. As a general feature, among the various time constants 
contributing to the damping, the more relevant ones correspond to the eigenstates that at 
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the same time possess larger energy separation and appreciable contribution to the state 
itself. 

Also here the WKB approximation allows one to estimate the damping time constants. 
By considering a rough picture in which a state is only made of eigenstates centered around 
n with width 2 Aw, all contributing with the same weight, the dynamics is ruled by the 
farest energy eigenvalues: 

2 



K E \E, 



n+An 
4/3 < 



— E, 



n—An 



2>Tl J KfiTl 



n + An 



1*2/3 



n — An 



1 \ 2/3 l 



(31) 



IV. PHENOMENOLOGICAL CONSIDERATIONS 

In this section we analyze observable consequence of the measurement of energy in the 
two cases of harmonic oscillators and bouncing particles in the presence of gravitational 
fields. Before doing this, some preliminar considerations are needed about the relevant 
coupling constant of the model, Ke- Bearing in mind that this parameter, describing the 
effective coupling of the apparatus to the observed system, intrinsically depends upon the 
particular experimental setup, one can only assess for it a reasonable value. If the analysis 
for a continuous nonselective measurement of energy is applied to the results of Itano et al. 
as done in |12|], a lower bound for the coupling parameter ke is found. In the coupling 



parameter is expressed in terms of a critical value K cr u = 4u;_r, where ujr = 12.272.s~ 1 is 
the Rabi angular frequency (corresponding to a period of the radiofrequency of 256 ms) and 
the energy difference between the two levels Ei — E\ = 2.125 ■ 10~ 25 J. The data of the 
experiment are fitted for K e ^ p w 10 2 K crit /(E2 — Ei) 2 = 10 53 J~ 2 s -1 . This value will be 
assumed, in some of the following examples, as an indicative one. 

It is worth to note that, due to the generality of the model, observable effects are in 
principle expected in all physical situations where quantum measurements are involved. 
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Among all the possible systems in which the effect can be made observable, we have chosen 
to discuss atomic or molecular systems and single degrees of freedom of macroscopic bodies. 

Let us first reconsider the case of systems modelizable as harmonic oscillators. If, for 
example, the vibrational energy levels of a biatomic molecule are monitored and the posi- 
tion operator Q(t) is interpreted as describing the istantaneous electric dipole moment of 
the molecule along the internuclear axis, then, according to flI?D, one can roughly expects 
a modification of the law according to which emission and absorption of electromagnetic 
radiation occurs. For instance, the intensity I(t) will not decay with an exponential law 
ruled by a lifetime r, but it will manifest a more complex behaviour: 

I(t) = Jo ex P { — ex P [^E{f iUJ ) 2 t]~ } (32) 

corresponding to an inhibition of the decay. Unfortunately, the time scale resulting when the 
K e g P quoted before is of the order of 10 8 s for reasonable vibrational frequencies, prohibitively 
long compared to the typical lifetimes of the vibrational transitions. Other perspectives can 
be opened by considering the time and frequency resolved spectra of spontaneous emission 



recently demonstrated [18 



Alternatively one could measure the energy in mechanical harmonic oscillators such as the 
resonators used as gravitational wave antennas, provided that the quantum limit is achieved 
in such a class of detectors [B5|,E0 . Due to the presence of a single system the analysis 



should be carried out by using the selective measurements approach. In this case, however, 
one big problem is the difficulty to obtain electromechanical transducers that measure the 
energy of the oscillator with enough sensitivity ||21|| . 

Other possibilities are also open by exploiting the quantum measurement model applied 
to the bouncing particle of Section III. An experiment aimed at testing the quantum damping 
can be designed on the basis of already performed experiments as the ones described in || 
and in 0. A cloud of atoms is trapped and cooled at a given height over a dielectric surface 
as in . Along the vertical path we put both a inhomogeneous magnetic field with constant 
gradient and a radiofrequency. The magnetic field is such that at the initial height the 
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radiofrequency gives rise to resonant Rabi oscillations between two levels of the atoms. The 
hyperfine splitting changes with the magnetic field and for each height less than the initial 
one the resonant condition is not fulfilled. A continuously operating laser acting along the 
vertical direction is tuned to the optical transition between level 1 and a third level which 
has forbidden transitions with level 2. A set of optical detectors allows one to observe the 
fluorescence light proportional to the occupancy of level 1. On the other hand, as shown in 



12], the occupancy of a level can be also thought as a measurement of energy since the state 
projectors of the occupancy and the energy operators coincide apart from a dimensional 
constant. Due to the spatially variable magnetic field, the continuous measurement of the 
occupancy will push the state of the atoms toward gravitational energy eigenstates and this 
will affect their average position. The damping of the average position of the atomic cloud 
could be measured with a destructive probe photon beam by repeating the measurements 
many times, as in |7|], or by exploiting non-destructive measurement schemes, such as the 



quantum nondemolition dispersive atomic probe one [^2J, allowing to repeatedly monitor the 
same atomic cloud (for more recent proposal see also f23| , |2"4|| ) . An estimate of the quantum 
damping time can be given by using Eq. (|3lD once expressed in terms of the more accessible 
position variance of the atomic cloud. The energy spread due to the z-motion of the atomic 
cloud can be written as 

AE 2 = (H 2 ) - (H) 2 , (33) 

where H = p 2 /2m + V(z). By supposing a Gaussian initial atomic phase space distribu- 
tion both in coordinate and momentum, centered on the z-axis at an height ze above the 
reflecting surface, we have a Wigner function 

W(z,p z ) = — !— exp [-{z - z E ) 2 /2z 2 ] ■ exp[-^/2 P y (34) 



and the evaluation of the energy spread gives in this case |16|] 



AE = mgz [2N\^y + (^-) 2 } 1/2 (35) 
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where the uncertainty product of coordinate and momentum in units of % has been intro- 
duced Zip z Jh = N. The energy spread has a value equal to AE ps mgzi for large values of 
Zi/zo, as one expects from the classical behaviour. In the opposite situation the effect of the 
Heisenberg principle appears and gives an inverse law dependence of the energy width upon 
the initial position spreading and a branching of the curves for different values of the un- 
certainty product N. A minimum value of the energy width is obtained for an intermediate 
value equal to z { = 2 1 / 3 N 2 / 3 z . 

The quantum decay constant is expressed as 

T= Z = ^-^[ 2N >(*L) i + (*)Y ( 36 ) 
k e AE 2 K E m 2 g 2 ZQ l K z { ' K z Q J J 

and in correspondence to the minimum of the energy spread gets the maximum value equal 
to r max (iV) = 2- 5 / 3 ^ 1 m"VV 2 iV- 4/3 . 

In Fig. 4 the decay constant versus the ratio z%/ zq is shown for various values of the 
normalized uncertainty product N, including the case of a pure state (N = 1/2). Smaller 
decay times are observed either for large values of the initial position uncertainty, in the 
right part of all the curves, or in the branching of the curves for the left part, in this 
last case depending upon the normalized uncertainty product. In the estimate we have 
assumed a constant width of the atomic cloud, an hypothesis that is not strictly valid due 
to the spreading following its preparation. To minimize this spreading a very low starting 
temperature of the atomic cloud is required, for instance 85 Rb clouds were shown to double 
their diameter in 15 ms if cooled at a temperature of about 10/iK [|25|] . Since low temperatures 
are also associated to small initial diameters of the cloud we expect the approximation of 
constant diameter to better hold in the left part of the curves drawn in Fig. 4. In these cases 
values of the quantum damping comparable to the ones of the right part are obtained for 
widths around ten-one hundred times the fundamental length zq, i.e. around 2 -=-20/zm. For 
comparison, the point corresponding to the decay observed in [J| is reported. To attribute 
the damping to the predicted quantum effect the classical sources of damping should be kept 
as small as possible, due also to the dynamics imposed by the Rabi transition frequency. 
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Currently achieved values of the escape time of atomic clouds are in the range of 10 seconds 
||26l , |27|1 and if this figure can be maintained together with a measurement coupling of the order 
of magnitude of the one corresponding to the experiment described in || and a normalized 
uncertainty product ~ 2 • 10 5 , quantum damping is made observable. 

V. CONCLUSIONS 

A quantum damping without classical counterpart has been introduced as a consequence 
of quantum measurements of energy and discussed for two situations. The model of the 
harmonic oscillator could be implemented in the monitoring of the vibrational motion of 
biatomic molecules or in quantum limited measurements in macroscopic mechanical res- 
onators. Quantum localization can also be manifested using a cloud of atoms bouncing over 
a reflecting surface in the presence of a uniform gravitational field. A possible experimen- 
tal scheme based on this last configuration has been discussed in more detail leading to a 
proposal that should merge two experiments already separately performed in || and 0. 
Studies of the quantized structure of particles in gravitational fields and the observation of 
a damping purely connected to the effect of the measurement as described here give further 
motivations to improve the cooling capabilities of atomic traps. 
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FIGURES 

FIG. 1. Average position versus time in the case of a particle bouncing on a reflecting surface 
for an initial pure state superposition of the first two energy eigenstates with amplitude coefficients 
1/2 and respectively. The case a) is relative to an unmeasured system (ke = 0) and the 

dashed curves represent the constant values for the two eigenstates <j>\ (below) and fa (above). 
The cases b) and c) are relative to Kg = 1CT 2 and ke = 10 -1 representing two examples in the 
underdamped and overdamped regimes, respectively. Here and in the analogous cases of Figs. 2 
and 3 we put H = m = 1. 

FIG. 2. Average position versus time for a particle bouncing on a wall and schematized by a 
Gaussian state whose center is initially located at h = 10 with variance a = 1. In the energy 
eigenstates expansion the main contributions stem from the eigenvalues between the fifth and the 
nineth. In a) the unmeasured case is depicted, in b) a measured case with ke = 10~ 2 is shown. Note 
the persistence of the oscillations which increase after a minimum indicating a not pure exponential 
damping unlike the one of Fig. 1. The lack of complete periodicity of the unmeasured case here 
and in the following Fig. 3a is attributable to the presence of many eigenstates contributing to the 
wavefunction reconstruction on timescales longer than the one depicted. 

FIG. 3. The same as in Fig. 2 but for a Gaussian state with variance a = 3. The larger 
spreading corresponds to an increase of the number of eigenstates which significantly contribute 
to the state, in this case between n = 3 and n = 12. Case a) is the unmeasured case, case b) is 
relative to a measurement with ke = 10 -2 . 
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FIG. 4. Decay constant for the quantum damping versus the position spreading corresponding 
to a Gaussian Wigner distribution of Cesium atoms bouncing in a gravitational cavity, normal- 
ized to the fundamental gravitational length zq (mass mcs = 2 • 10 _25 Kg, gravitational length 
zq = 0.23//m). The curves are obtained for different values of the normalized uncertainty product 
of the atomic cloud, respectively N = 1/2 (pure state, a), N = 20 (b), N = 2 • 10 3 (as experimen- 



tally achieved in @, c), and N = 2 ■ 10 5 (d). For comparison the experimental point explored in 
0], corresponding to a decay time of ~ 80 ms for a value of Zi/zQ ~ 10 3 and explained in terms 
of classical sources of damping, is shown. It has been assumed a measurement coupling constant 

K E = 10 53 J- 2 S _1 . 
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